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Abstract : 

The universal field equations introduced by the author and his collaborators, which 
admit infinitely many inequivalent Lagrangian formulations are shown to arise as consis- 
tency conditions for the existence of non-trivial solutions to the quasi-linear equations, 
called equations of hydrodynamic type by Novikov , Dubrovin and others. The solutions 
in closed form are only implicit. A method due to Stokes, which is in essence just Fourier 
Analysis is resurrected for application to those equations. With the benefit of algebraic 
computation facilities, this method, allows the general structure of power series solutions 
to be conjectured. 
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1. Introduction. 

The purpose of this article is to comment upon some features of nonhnear equations 
related to the Euler equations of motion for fluid flow, which have attracted consider- 
able attention among theoretical physicists in recent months after the seminal paper of 
Polyakovl^l Two diverse lines of inquiry regarding the properties of these equations are 
followed; The flrst is to demonstrate how the equations which we have called Universal 
Field Equations [^~^] may be derived as consistency conditions on the integrability of a 
simple set of first order equations which have been extensively investigated by Dubrovin 
and Novikovi^'^] and by Tsarev^^"^^! under the name of 'systems of hydrodynamic type'. 
The Universal Field Equations enjoy remarkable properties. They may be derived from 
an infinite number of inequivalent Lagrangians, are associated with an infinite number of 
conservation laws and are either reparametrization invariant t^^, or else the solutions are 
completely covariant. They may be linearised by a Legendre Transform and thus integrated 
to give an implicit solution!^]. 

The second is concerned with an iterative method of solving such nonlinear equations 
which goes back to Stokes t-*^-*^], as referenced in Whitham's bookl-*^^!. Stokes did not have the 
advantage of a sophisticated algebraic manipulation package, such as MAPLEt-*^^! which 
permits the rapid evaluation of approximations to a sufficient order to admit their structure 
to be guessed, otherwise the development of methods of dealing with integrable equations 
might have taken a different course. Consider first the case of (1+1) dimensions; take 
the simple nonlinear wave equation, ut = uux, or perhaps better the KdV equation ut = 
6uUx + Uxxxj OT the Burgers equation, ut = uUx + i^Uxx, which is the Navier Stokes equation 
in the absence of a pressure gradient in one space dimension. All of these equations can be 
solved in principle; the first by the method of characteristics to give the implicit solution 
u — F{x + ut), the second by the Inverse Scatttering method and the third by the Cole- 
Hopf transformation. In practice an explicit approximate solution might be more useful, 
as frequently the problem posed is that of the time evolution of an initial wave-form. Look 
for a solution of the form of a formal Fourier expansion 

u{x, t) = uo{t) + ui{t)e''^ + W2 We'"" + i^s^e'"" + . . . . (1.1) 

The initial data is coded into the values of the coefficients Uj{t) at t = 0. Then the key 
point is that these coefficients are determined recursively by a set of inhomogeneous linear 
equations, where the driving term for the equation for Un{t) is determined by terms of the 
coefficients Uj{t), j < n. uq is a constant and Ui is determined by a homogeneous linear 
equation; U2 is then determined by an inhomogeneous linear equation with driving term 
determined by ui and so on. The general structure of the solution takes a remarkably sim- 
ple form, and allows a Taylor like expansion for solutions to the non-linear wave equation 
to be formally deduced. If one looks for periodic solutions of the KdV equation, then this 
method gives a sequence of approximations to a given initial value problem. This is just 
what the inverse scattering method was set up to achieve! Admittedly this construction is 
formal; but then the inverse scattering method permits a exact solution in closed form to 
the initial value problem, only when this corresponds to a collection of of solitons! Indeed 
the multisoliton solution of the KdV equationt^^l can be recovered by formal summations 
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using the above methodf-'^^l . Another possible apphcation of this approach is to the quanti- 
sation of such nonhnear equations, by considering the coefficients in the expansion of (1.1) 
as creation operators. 

2. The Universal Field Equations as Consistency Conditions. 

Consider the following set of equations 



AT 



dUi{Xk) n • 1 AT 



i=i 



(2.1) 



A special case, in which u{l) = — 1, Xi = t constitutes a multidimensional version of the the 
equations of hydrodynamic type^^^ ,or equivalently, (2.1) may be regarded as a static limit 
of those multidimensional equations. For the purposes of this section it is more convenient 
to treat all the variables on the same footing, as in (2.1). Now suppose all Uj, j = 2, . . . , N 
are all functions of a single function ui = v. Then the result of differentiating (2.1) for 
i = 1 with respect to xi, 3:2, . . . xn in turn and using the easily proved result that 
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yields the N equations 

UlVx^X2 + U2VX2X2 + 



+ UNVxiXN+i 

+ UNVx2XN + i 



dm 



'^I'^xiXn + U2VX2XN + ■ • ■ + UnVxj^xn + { 



dxi 


+ 


dx2 


dui 




du2 


+ 




dxi 


dx2 




: 




dui 




du2 




+ 





dx\ dx2 



+ ... + ^ — )^;i = 

OXn 

+ • • . + ^ )V2 = 

OXn 

+ • ■ • + = 



(2.3) 



dx 



N 



Elimination of the + 1 variables Ui, i = 1, . . . N and between the N + 1 equations 

(2.1) and (2.3) yields the determinantal condition 



det 



— 

dxk 



dv 



= 



(2.4) 



dxjdx^; 



This is just the Universal Field equation of reference [1]. Furthermore, the extension to 
incorporate many fields arises in a similar manner. Suppose the only the first n among the 
N functions Uj are functionally independent; i.e. the functions Uj, j — n + 1, . . . , N are 
functions of Ui, i = 1, . . .n. Then consider the equation 



N 



fc=0 



dx-i 



(2.5) 



3 



Here A^, k = 0, . . . ,n are a set of arbitrary constants. Differentiation with respect to Xj 
gives 




Sr^ v-^ d Uk , \^ .sr^ . duk , -r-^ oUr sr^ . duk. „ 

2:-.E^^a^ + 2:a^E^'=a^+ E a^2:^^a^]=o- (2-6) 

j=Q k=0 ■' j=0 k ■' r=n+l ■' k=0 

Elimination of the N + n variables Uj and 

(Efe Afc|^+Ef=n+i 1^ ELo ^fefe) from the N equations (2.6) regarded as linear equa- 
tions in those variables and the n independent relations among the equations 

reproduces the multifield form of the Universal equation. In compact form, these gener- 
alised Bateman equations are given by 

dUg 

= . (2.8) 

Here, the determinant is that of a (N + n) x (N + n) matrix, with n being the number of 
fields Ua and N the number of coordinates Xi {n < N), and the indices (a, i) and {b,j) refer 
to rows and columns respectively (in particular, the entry "0" stands for the n x n null 
matrix). The quantities A'^ are arbitrary coefficients, in terms of which the determinant is 
to be expanded. It is understood that (2.8) has to hold for all values of these coefficients. 
Hence, one obtains (^Zi) equations - clearly generalising the Bateman equation (2.4) 
for one field to many - whose general covariance under arbitrary field redefinitions in u°' 
is easily established. Even though these equations form an over determined set, except 
for (n = 1) or (n = N — 1) (the equations in the latter case were shownt^l to be also 
universal), their space of solutions is non empty I^^, as any arbitrary choice of functions 
Ui{xj), homogeneous of degree zero in the variables Xi will satisfy (2.8). 

3. Introduction of explicit time. 

Consider the equations, 

ut + uux + vuy =0; 

Vt + UVx + VVy =0. 

These equations are related to the one Polyakov and others have studied in the theory 
of conformal turbulence through setting u = i)y\v — —ipx- This constraint ensures that 
the two dimensional velocity vector {u, v} describes incompressible fiow; i.e. Ux + Uy = 0. 
Of course the above equations simply say that the total derivative ^u{x, y, t) — 0, and are 
the Euler equations for the free motion of an inviscid fluid. The equation for the vorticity 
a;, which is given hy uj — Uy — Vx — V'^i^ which Polyakov studies is 

LO = IpxU^y - i>yOJx (3.2) 
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There is a further consistency condition on 7/;, arising from (3.1), namely 



V'xx^yy - Iply = 0. (3.3) 

if in addition to (3.2) equations (3.1) are imposed. These considerations have been extended 
to the Navier Stokes equationl^^l by the intoduction of a viscous force in (3.2) Somewhat 
curiously, if an irrotational flow with a; = is sought, then {u, v} is derivable from a 
potential (j) as V0, and the equations (3.1) reduce to the following single equation for 0, 

0t + ^W)' = O. (3.4) 

Of course, the requirement of incompressibility imposes an additional constraint upon 0, 
namely that it should satisfy the equation = 0. It is easy to see, by transforming to 
complex variables z — x + iy, z — x — iy, that the only solution in this case is {u, v} = 
constant. There are infinitely many Lagrangians; 

f ( UV] 

ut + f{uv)ux - f{uv)vy. (3.5) 

An implicit solution; 

u =F(x — ut,y — vt) 
V =G{x — ut,y — vt). 

where F and G are arbitrary functions of two variables can easily be verified, and this 
result extends in an obvious manner to the solution of the equations 

dui dui , 

That these equations admit solutions which develop singularities can be seen by looking 
for solutions which satisfy the ansatz 

uAt,x,) = ^^ (3.8) 

With the choice Zi{xk) = J2 -^ikXk, then (3.8) is a a solution provided the matrix A satisfies 
A = A'^A. 

Introduce another set of fields Vi, i = 1 . . . A^.Then 

j=i j=i fc=i 

is a Lagrangian for the equations (2.4). The companion set of equations, obtained by 
varying jC with respect to Ui is 

dvi J^ . dvi dun duj . 



If u is eliminated from equations (3.1) the resulting equation for v can be written in 
a nice determinant al form: 
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(3.11) 



Note that if v docs not depend upon y this equation is the Bateman equation, various 
aspects of which are considered in references [2-5] .In that case the complete solution is 
given by the implicit solution for v of the equation 



tf{v) + xg{v) = c 



(3.12) 



where /, g are arbitrary functions and c is a constant which may be zero. In the general 
situation the solution is given by the elimination of u from (3.6) 



4. Explicit Formal Solution of the Nonlinear Wave Equation. 

It is well known that the general solution of the nonlinear wave equation 



Ut = UUa 



is solved implicitly as 



u = G{x + ut). 



[12] 



(4.1) 

(4.2) 



It may however be solved for an arbitrary differentiable initial value F in terms of a Laplace, 
or Fourier expansion of the initial value of the form 



F{x) = Co + cie"-"" + C2e^ax + cse^"^ + 



(4.3) 



by the method outlined in the introduction, cq may be taken as 1, by rescaling t. A com- 
puter calculation using MAFLEt-*^^] reveals the conjecture that the power series expansion 
of u{x, t) is given by 



u 



sr^ '^c^"^(Ari)"-iexpiV(ax + i) 

ni,n2,... j=\ 



(4.4) 



Here N — rii + 2n2 + ■ ■ ■ kuk, n — ni + n2 + ■ ■ .Uk- These coefficients may displayed as a 
formal contour integral, and the formal sum constructed to give 
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(4.5) 



=F(0) + 5^f 



n=l 



dx'^- 



after rescaling so that F{0) = cq is arbitrary. The contour is taken sufficiently large so as 
to encircle the multiple poles ai z = t The sum may itself be formally performed, to give 
the result 



1 /, tY.^CjCxi)aji.r + z)-^dz 



=F{0) - ^ /log(l - J^^iFix + F{0)z) - F(0))^ 



(4.6) 



There does not appear to be a simple proof of the result (4.5), but it may be used to obtain 
the explicit solution 

u{x,t) = ^, (4.7) 

which also follows fom (4.2). It would be interesting to extend this to (3.1). We may note 
that if v{t) is any prescribed function of u{t) the two equations (3.1) reduce to a single 
equation, 

ut + uux + f{u)uy = 0, (4.8) 

where / is arbitrary. This case at least, is not tractable by the above method. However, 
the physically interesting case is when u, v are derivable from a stream function, and 
Ux + Vy = 0. An explicit, possibly general, easily verifiable class of solutions dependent 
upon an arbitrary function of a linear combination of x, y, t can be found by the Fourier 
method as 

it(x, y, t) =uo + bG(ax + by+ (auo + bvo)t), 
v{x, y, t) =vo — aG{ax -\-by + {auQ + bvo)t). 



5. Conclusion 

Some properties of the quasi linear system of equations of hydrodynamic type have 
been explored, and their general solution given in terms of arbitrary functions, but only in 
implicit form. This is a feature of some nonlinear completely integrable systems in higher 
dimension, and results from the fact that one of the methods of solution is to interchange 
the role of dependent and independent variables, whereupon the equation linearizes!^! How- 
ever, if a Fourier series solution, with time dependent coefficients is sought, then instead of 
the complete decoupling of the ordinary differential equations describing the time evolution 
of the modes, these are determined successively in terms of the fundamental mode. It is a 
remarkable feature that the structure of the power series solution to the simplest nonlinear 
wave equation is easily elucidated with the aid of an algebraic computation package, and a 
formal power series solution obtained for an arbitrary initial configuration for the simplest 
nonlinear wave equation. 
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